Origin of the Growing Length Scale in M-p-Spin Glass Models 
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Two versions of the M-p-spin glass model have been studied with the Migdal-Kadanoff renor- 
malization group approximation. The model with p = 3 and M = 3 has at mean-field level the 
ideal glass transition at the Kauzmann temperature and at lower temperatures still the Gardner 
transition to a state like that of an Ising spin glass in a field. The model with p = 3 and M = 2 has 
only the Gardner transition. In the dimensions studied, d = 2, 3 and 4, both models behave almost 
identically, indicating that the growing correlation length as the temperature is reduced in these 
models - the analogue of the point-to-set length scale - is not due to the mechanism postulated in 
the random first order transition theory of glasses, but is more like that expected on the analogy of 
glasses to the Ising spin glass in a field. 
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One of the leading contenders for a theory of glasses is 
the random first-order transition theory (RFOT) [TH3]. 
It had its genesis in p-spin glass models pQ . The particu- 
lar p-spin models which might be relevant to the proper- 
ties of structural glasses have a mean-field limit in which 
there are two critical temperatures Td and Tk- The up- 
per temperature Td marks the temperature at which dy- 
namical singularities appear and are like those found in 
mode-coupling theory [4]. The lower temperature Tk is 
the temperature of the ideal glass transition. This occurs 
where the configurational entropy (or log of the number 
of metastable states) vanishes [5] . Mean-field like calcula- 
tions on glass-forming liquid models support this picture 

mm- 

It has always been recognized that the dynamical tran- 
sition at Td will disappear outside the mean-field limit 
due to activated processes out of the metastable states. 
These activated processes make even the existence of 
metastable states problematical. In a recent paper Franz 
et al. [8] calculated a dynamical correlation length us- 
ing a field theoretic approach and found on comparing 
with numerical data that the agreement was only good 
when the length scale was of the order of a particle di- 
ameter. The field theory predicts that this length scale 
should diverge but the simulations reveal that the corre- 
lation length remains small [9], even though time scales 
increase rapidly: the dynamical transition is an example 
of an avoided transition. 

The ideal glass transition at Tk at mean-field level (or 
infinite dimension) is associated with a static (equilib- 
rium) transition to a state with one-step replica symme- 
try breaking (1RSB) pQ. The order parameter q jumps 
from zero in the high temperature phase to a finite value 
at and below Tk- It is this jump in q which leads to the 
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"first-order" part in the name of the RFOT theory. The 
configurational entropy (or complexity) vanishes as the 
temperature approaches Tk from above. While there is 
widespread agreement that the transition at Td becomes 
just a crossover or avoided transition in finite dimensions, 
there is no consensus about what happens to the ideal 
glass transition outside the mean- field limit. One of us 
has argued that the 1RSB transition must also be avoided 
in any finite dimension [lOj. just like the dynamical tran- 
sition at Td- In other words, the lower critical dimension 
d of the 1RSB state is infinite. 

In this paper we examine M-p-spin models; in partic- 
ular the cases of p = 3 with M = 2 and M = 3. These 
variants of the p-spin model have been extensively stud- 
ied [THUS] . Their significance is that calculations and 
simulations can be done with them both at the mean-field 
level and in finite dimensions. In the M-p-spin model, 
there are M Ising spins a^ a \ a = 1, 2, • • • , M on each 
site i of (say) a hypercubic lattice. The spins interact 
with each other via a p-body interaction. The Hamilto- 
nian involves terms of products of p spins chosen from the 
spins in a pair of nearest-neighbor sites. For the p = 3 
case, the Hamiltonian is given by 

M M 
(ij) a</3 7 

+4^a^af), (1) 

where the notation (ij) means that the sum is over all 
nearest neighbor pairs i and j. The number of different 
coupling constants, j}"^ 7 and Jj^ 01 ^ for given (ij) is 

just 2M(^f) = M 2 (M - 1). All these couplings are usu- 
ally chosen independently from a Gaussian distribution 
with zero mean and width J. 

The versions with M — 2 and M — 3 are of particular 
interest as at mean-field level they have quite different 
kinds of behavior. The model with M — 3 at mean- 
field level has both a dynamical transition at Td and an 
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ideal glass transition at Tk- The model with M = 2 is 
completely different at mean-field level. It has neither 
of these transitions. The origin of the differences can 
be glimpsed by putting the M-p-spm model into a field 
theoretical framework. 

The standard way of doing this is to use the Hubbard- 
Stratonovich transformation on the replicated partition 
function and then trace over the spins. The resulting 
field theory associated with this model is the following 
Ginzburg-Landau- Wilson Hamiltonian 

n GLW = /d d r[^[Vg a6 (r)] 2 +*^^(r) 

k a<b a<b 

a<b ) 

where <? a &(r) is the order parameter and a and b are 
replica indices running from 1 to n with n — > 0. At mean 
field level, this model has been known for a long time to 
show very different behavior depending on the value of 
R = W2/1V1 [IB]. When R > 1, there are two transitions 
at the mean field level as described above; a dynamical 
transition at some temperature and a thermodynamic 
transition at a lower temperature Tk to a state with one- 
step replica symmetry breaking. When R < 1 neither of 
these transitions will occur. In Ref. [17], the ratio R 
was evaluated for the M-p-spin model for general values 
of M and p. The cases we are interested in this paper, 
namely p = 3, M = 2 and p = 3,M = 3, correspond 
to R m 0.879 and R — 2, respectively. Therefore the 
two models should indeed show very different mean-field 
behavior. 

At temperatures below Tk for the case where R > 1, 
there is yet another transition, the Gardner transition, to 
a state with full replica symmetry breaking (RSB) [18]. 
For the case where R < 1, there is at mean- field level 
only one transition - the Gardner transition - again to 
the state with full RSB. 

The transition discovered by Gardner is thus present 
for both the M = 2 and M = 3 models at mean-field 
level [18]. She showed that the state with full replica 
symmetry breaking was very similar to that of the Ising 
spin glass in an applied field h. For this model, there is a 
line in the h — T phase diagram, the de Almeida-Thouless 
(AT) line, [19], which separates the paramagnetic replica 
symmetric state from the state with full replica symmetry 
breaking. Arguments have been presented [SBJ [UJ that 
the lower critical dimension of states with full replica 
symmetry breaking is 6. The Gardner transition, which 
is in the same universality class as the AT transition, 
should be another avoided transition for all d < 6. 

In this paper we have studied both the models with 
M = 2 and M = 3 within the Migdal-Kadanoff (MK) 
renormalization group approximation in dimensions d = 
2, 3 and 4 to determine how thermal fluctuations modify 
the mean- field picture of these two models. The MK 
approximation is one of the few approximations which is 



TABLE I: The correlation length £(0) in the zero-temperature 
limit for various values of M and d measured in units of the 
lattice spacing. 
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reliable for the study of spin glasses in low dimensions 
[IT] [T2] . The details of our calculation are as in [IT] [T2] . 
We are interested in particular whether in the physically 
relevant dimensions, d — 2 and d = 3, whether there 
are any vestiges left of the mean-field transitions. One 
can see in the molecular dynamics study of Kob et al. 
[3], clear remnants of the dynamical transition. Only 
equilibrium properties are studied in this paper so only 
the remnants of transitions which could be seen are those 
of the ideal glass transition and the Gardner transition 
for the case with M = 3, and just the Gardner transition 
for the case M = 2. 

We determined the correlation length £ by the same 
method which was used in Refs. [UJ [T2l [15], that is 
via the decay of the interactions with distance L on 
the MK hierarchical lattice: Jij ~ exp(— L/£). As the 
temperature T is reduced to zero this correlation length 
grows to a value £(0), which is strikingly large especially 
for d = 4. The data on £(0) are presented in Table]!] The 
large values of £ (0) certainly suggests there is an avoided 
transition mechanism at work. Even in d = 3, £(0) is 
considerably larger than those which have been obtained 
in simulations of realistic glass models, at least down to 
temperatures which are currently practical [9 . 

It is useful to measure temperature T on the scale 
of the mean-field transition temperature T c (defined as 
when the coefficient t in Eq. ^ equals zero). For M = 2, 
T c = (V2z)^ 2 J, while for M = 3, T c = (3z)^ 2 J, where 
z = 2d is the number of nearest neighbors on the hyper- 
cubic lattice [T7]. In Fig. [lithe ratio of £/£(0) has been 
plotted as a function of T~jT c . It shows that as a func- 
tion of T/T c , the ratio £/£(0) is essentially the same for 
both M = 2 and M = 3. We had expected to see for the 
case of M = 3 features which could be associated with 
a possibly avoided ideal glass transition at Tk- None is 
visible in Fig. [I] This result is our main finding. What 
is its significance? 

The correlation length studied here is the equivalent of 
the point-to-set length scale [HI [55] [53] studied in glassy 
supercooled liquids. In the RFOT theory [H[5] this length 
scale grows as the temperature is reduced and eventually 
diverges at Tk- In that theory the growth of the corre- 
lation length is driven by the decrease of the configura- 
tional entropy (or complexity) to zero as the tempera- 
ture approaches Tk- Since for M = 2 there is no ideal 
glass transition yet this model is almost identical in its 
properties with those of the M = 3 model, the growing 
correlation length cannot be a remnant of the ideal glass 
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FIG. 1: Behavior of the scaled correlation function £/£(0) 
plotted versus the scaled temperature T/T c . The closed sym- 
bols correspond to M — 3, the open symbols to M = 2. The 
upper (green triangles) curves are for d = 4, the middle curves 
(blue squares) are for d — 3, while the lower two curves (red 
circles) are for d — 2. 



transition. It must instead be a remnant of the Gardner 
transition which is common to both models. 

Thus the mechanism behind the growing correlation 
length as the temperature is reduced cannot be that en- 
visaged in the RFOT, but instead must be that associ- 
ated with the growing correlation length which arises in 
the Ising spin glass in a field as the temperature is low- 
ered. According to the droplet picture, .24-26] the cor- 
relation length increases as the temperature is decreased 
and saturates at T = to a value set by equating the 
interface energy between a droplet of size £(0) and its 
time reverse, ~ J£(0) e , to the energy gained from the 
field on flipping the droplet, ~ ft,£(0) d / 2 . The exponent 
9 « -0.28 for d = 2 and 9 w 0.24 for d = 3 (see Ref. 
[27j for a review of the value of 9 in various dimensions 
d). Table 1 shows that £(0) gets larger as the dimension- 
ality goes up and when there is an AT line, i.e. when 
d > 6, it would be expected to be infinite. The MK ap- 
proximation itself is a low-dimensional approximation (it 
is exact in one dimension), and cannot be trusted to be 
even qualitatively correct in dimensions as high as six. 

We have argued before that the growing (point-to-set) 
length scale seen in supercooled liquids as the tempera- 
ture is reduced [351 US] is a consequence of their being in 
the same "universality class" as the Ising spin glass in a 
field. However, until now, we could not rule out the pos- 
sibility that the growing length scale might arise through 
a 1RSB transition as in RFOT, (but possibly avoided 
because of the mechanism in Ref. [ID]). The similarity 
between the M = 2 and M = 3 models shown in Fig. [I] 
now removes that possibility. 

In RFOT theory, the configurational entropy is argued 
to go to zero at the Kauzmann temperature Tk- For 
M-p spin glass models it is not clear how the configura- 
tional entropy should be determined outside the mean- 



FIG. 2: The entropy S per site as a function of the scaled 
temperature T/T c . The solid lines are for the M — 3 model, 
the dashed lines are for the M — 2 model. 
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FIG. 3: The Edwards- Anderson order parameter q as a func- 
tion of the reduced temperature. The symbols have the same 
association with M and d as in Fig. [l] 

field limit, but we have studied their total entropy by nu- 
merically differentiating the free energy calculated within 
the MK approximation (which leads to some inaccuracy 
near T = 0). It is plotted in Fig. § as a function of 
the reduced temperature T/T c . Once again the models 
with M — 2 and M = 3 behave almost identically, indi- 
cating that when the correlation length gets large, there 
is present a form of universality. At high temperatures 
where the correlation length is small, the two types of 
model have very different entropies: the high tempera- 
ture limit of the entropy per site is feM In 2. While there 
is no sign in Fig. [2] of the entropy vanishing below some 
temperature Tk , the entropy is smaller at lower tempera- 
tures for the d = 4 versions of the model, and behavior in 
four dimensions is going to be closer to mean-field theory 
than behavior in two dimensions. 

We have also determined the Edwards- Anderson order 
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parameter q 

q = [K (Q) > 2 U (3) 

where the average is over the bond realisations, q is in- 
dependent of the value of a, (which runs from a = 1 
to M). Under the MK iteration scheme the couplings 
flow to the high-temperature fixed point where the block 
spins are decoupled and only single site terms remain. 
As a consequence it is easy to evaluate q. 

The Edwards- Anderson order parameter q is plotted as 
a function of the reduced temperature in Fig. [3j The fig- 
ure shows that q —> 1 as T — >• which is to be expected. 
At low temperatures where the correlation length is large 
both the M = 2 and M = 3 models behave almost iden- 
tically, which is another example of the "universality" 
emerging in the problem. 

What is striking is that q is non-zero at any tempera- 
ture, although it does become very small when T>T C . 
q is a measure of the extent to which the spin at site i 
remembers its initial orientation: i.e. 

1=^H(^ a) (^ a \t)), t^co. (4) 

i 

Thus in the p-spin models the spins never completely 
forget their initial orientations, no matter how high the 



temperature. This behavior is not a consequence of using 
the MK approximation. It is a feature which arises in 
any model described by the "Hqlw 01 Eq. p]) with a 
non-zero value of W2, the term which breaks the time- 
reversal symmetry. At mean-field level q does vanish at 
temperatures above T^. 

p-spin models are meant to be useful models for un- 
derstanding the properties of supercooled liquids so this 
feature of a non-vanishing q is hard to reconcile with 
the properties of supercooled liquids. These forget their 
initial conditions after a time of the order of the alpha 
relaxation time, so for them q is zero on long time scales. 
Maybe p-spin models are useful for describing the prop- 
erties of supercooled liquids but only on timescales less 
than the alpha relaxation time. Given the huge effort 
which has gone into investigating their properties, this is 
certainly to be hoped. 

Finally we note that if we had used the MK approxima- 
tion with the further approximations which were made in 
Ref . [30j , a Kauzmann transition would have been found 
for the model with M — 3. It is only by carrying out 
the MK calculation exactly that one recovers the correct 
behavior [TT| . 
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